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Abstract. Massive particles are radiated from black holes through the Hawking mechanism to- 
gether with the more familiar radiation of massless particles. For E >m, the emission rate is iden- 
tical to the massless case. But E < m particles can also tunnel across the horizon. A study of the 
dispersion relation and wave packet simulations show that their classical trajectory is similar to that 
of a boomerang. The tunneling formalism is used to calculate the probability for detecting such 
O , E < m particles, for a Schwarzschild black hole of astrophysical size or in an analogue gravity 

experiment, as a function of the distance from the horizon and the energy of the particle. 
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^ INTRODUCTION 

(N 

^ The standard vision on Hawking radiation (HR) of massive particles dates back to [1]. It 

• states essentially that there is a cut-off in the emission of massive particles because 

(^ the Hawking temperature Th of the black hole must deliver E > m. Therefore, as 

CN anticipated in [2], "there will not be much emission of particles of rest mass m unless the 

Vl temperature k/Itt is greater than m". The example of an electron gives nie ~ I0~^'^kg, 

^ >> hence T ~ 10^.^ and so the black hole mass must decrease to M ~ lO^^^niQ ~ Kfi^nipi 

S^ {niQ and nipi are the solar and Planck mass, respectively) before electrons can be emitted 

H in any significant amount. Through the emission of massless particles, the black hole will 

evaporate and this threshold eventually be reached. However, during most of the lifetime 
of an astrophysical black hole, there will be no significant emission of massive particles. 
This standard vision sees HR as a "black box" whose outcome is measured at asymp- 
totic infinity. It must be amended if one views HR as a near-horizon process. This allows, 
e.g., to reconcile the two heuristic interpretations of HR as mode conversion between 
positive/negative-norm partners, related to pair creation (from vacuum fluctuations) just 
inside the horizon, or just outside the horizon. Obviously, any calculation should give the 
same result in both interpretations. A thermal spectrum is indeed detected at asymptotic 
infinity in both cases. But the mass decrease of the black hole is due to the emission of 
the positive partner in the first case, and to the absorption of the negative partner (in the 
second). There might be good reasons for the emission of positive partners to be subject 
to the threshold E > m, but there is no reason why the absorption of negative partners 



should obey a similar condition. How can these two heuristic visions then be equivalent? 

DISPERSION RELATION 

The simplest dispersion relation for massive modes/particles is (in units h— I) 

{(0-Uk)^ = m^ + c^k^, (1) 

where t/(r) < is the free-fall velocity of an observer starting at rest at infinity. 



U{r) = —c\/rhlr for a Schwarzschild black hole (with r/, the horizon or Schwarzschild 
radius, i.e. lJ{r}^ = —c), and we have defined the "mass" m = mc^. Such a Klein-Gordon 
dispersion relation can be derived for example from the Painleve-GuUstrand-Lemaitre 
form of the Schwarzschild metric 

ds^ = [c^ - U{rf]dt^ - 2U{r)dtdr - dr^, 

which provides an intuitive analogy with sound or surface waves counter-propagating 
with local velocity c against a background fluid flowing at a speed U in the laboratory 
frame. (Oq = \Jm?- + fc^ is then the frequency in the co-moving reference frame, and 
CO = (Of) + k ■ U is the Doppler-shifted frequency (corresponding to the laboratory or 
"black hole rest frame"). Note that ft) is a conserved quantity (for stationary spacetimes), 
and therefore provides a good basis for a semiclassical treatment (see below). 

A graphical analysis of the mode conversion characteristic of Hawking radiation for 
this dispersion relation is given in Fig. 1 [3]. The key point here is that the appearance 
of the negative mode solution (negative co-moving frequency, i.e.: bottom part of the 
figure) is identical for the various upper parts of the figure, i.e.: the positive-negative 
mode conversion around the horizon is totally independent of the ratio E/m (there is no 
threshold E > m) and in fact identical to the case of a massless mode [3]. 

It is also easy to realize that there exists a critical value U* for the counterflow in 
the case E < m, such that there are two positive mode solutions when \U\ > \U*\, 
and there is none left when \U\ < \U*\. U* constitutes a turning point or saddle-node 

bifurcation and its value is U* = icy 1 — (^) [3]. This shows that the E < m modes 
cannot escape to infinity but eventually turn around like a boomerang at a location 
r* = rh{c/U*)^ (which depends on their ratio E/m through U*) and fall back into the 
black hole (see [3] for a wave packet simulation). Therefore, they do not contribute 
to the evaporation of the black hole. Nevertheless, they can (at least in principle) be 
detected at any finite distance from the horizon. Putting in numbers gives r* ^ Sr^ for 
E = 0.8m, or r* ^ SOr^ for E = 0.99m. Taking into account that, for a solar-mass black 
hole, r/, ~ 3km, and moreover that the resolution of any detection will be limited by the 
time until re-absorption in the black hole, it seems that this boomerang effect only has a 
purely theoretical interest for astrophysical black holes (except perhaps for supermassive 
ones). 
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FIGURE 1. Positive and negative mode solutions for the dispersion relation (1) for different values of 
the "counter-current velocity" U and (a) E <m; (b) E > m. The mechanism of appearance of negative- 
norm solutions (bottom) is independent of the ratio E /m, and actually equal to the massless case [3]. 



TUNNELING FORMALISM 



The tunneling probability V7 in a semiclassical formalism h W °^ exp[— 2Im5], 
where Im5 = Im /p^(r)(ir is the imaginary part of the action along the semi- 
classical trajectory, and Pr{r) is obtained from the dispersion relation (1), which 
we now write (£ — p-U)^ = m^ + p^. The expected final result is of the form 
W{E) oc exp[— 2Im5i]exp[— 2Im52], where 5i corresponds to tunneling through the 
horizon, and 52 stems from the second classically prohibited region beyond the horizon 
in the case E < m. Recall that for U{r -^ oo) = 0, E^ = m^ + p^ (i.e., the standard 
Minkowski dispersion relation), and hence particles with £ < m are certainly classically 
forbidden at flat asymptotic infinity. From the dispersion relation, p can be written 
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\U^-l)+E^ = pi+P2. 



Si is obtained from p\ in the usual way [4]. We shift the contour of integration to the 
complex plane, apply the residue theorem Jcf{z) = IniResf, find a pole along the 
radial path: U = —I at r = r^, and obtain Respi = Resi 
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'{l+U){l-U) ^ 2U'{rh)- ^^^^ ^^^^^ 

to the standard Hawking result: Wi{E) ©« exp[— E/T/f], with Tfj = |f/'(r/,)|/2;r, where 
the prime denotes d/dr. This confirms that the presence of a mass has absolutely no 
influence on the probability of tunneling across the horizon. 

We now consider a detector at a distance R and examine the second contri- 
bution 52. There will be an imaginary contribution if E^ < nP-{\ — U^) (i.e.. 
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FIGURE 2. Semiclassical barrier lmp2 (thick lines) and tunneling action Im52 (dashed lines) beyond 
the horizon r/, for various values of E/m < 1 [5]. 



1/2 



E < Ec{R) = m{\-'^) ' for given R, or R > rc{E) 
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E). Ffor i? — )■ 00, we recover the condition E < m. For a Schwarzschild pro- 
file, we obtain Im52 = /,. Jr ^_^ , ^rrP-{l—rh/r)—E^, see Fig. 2. In the limit 
R ^ r/,, the integral is dominated by the contributions [/ — ;■ and one can write 
Im52 = \/rrP- — E^ J^ dr k, Ry/rn^ —E^ (assuming R ^ re). The overall tunneling rate 

then becomes 
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Other simple analytical results can be obtained in the near-horizon limit R^>- r^ [5]. 

ANALOGUE GRAVITY 

As mentioned above, the emission oiE <m particles has little practical impact for astro- 
physical black holes. However, in the case of analogue-gravity systems with a massive 
dispersion relation, they become more interesting. Indeed, in such systems, particles 
can be detected inside the black hole also, so there is no limit on the spectral resolution. 
Moreover, the region where U ^0 can be extended far away from the horizon. Candidate 
systems include acoustic waves in ion rings, massive phonons from 2-component BECs, 
Langmuir waves in a moving plasma, barotropic waves (inertia-gravity or Poincare 
waves), and spin waves in magnetic media [3]. Note that the existence of two imagi- 
nary contributions for E < m indicates the presence of a double barrier, and hence the 
possibility of creating resonant states according to the Bohr-Sommerfeld quantization 
condition. These resonant states could be interesting candidates for detection and hence 
for a confirmation of some of the curious features of Hawking radiation described here. 
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